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1 Introduction 


The CP-violating phase (ps arises in the interference between the amplitudes of mesons 
decaying via b —)■ ccs transitions to CP eigenstates directly and those decaying after 
oscillation. The phase (ps can be measured using the decay P° —)■ J/xjjcj). Within the 
Standard Model (SM), and ignoring penguin contributions to the decay, is predicted to 
be —2/3^, with j3s = arg(—where Vij are elements of the CKM matrix j^. 
The phase 0s is a sensitive probe of dynamics beyond the SM (BSM) since it has a 
very small theoretical uncertainty and BSM processes can contribute to B^-B^ mixing 
[^. Global fits to experimental data, excluding the direct measurements of (ps-, give 
—2/3s = —0.0363 ± 0.0013 rad [^. The current world average value is 0s = —0.015 ± 
0.035 rad Q, dominated by the LHCb measurement reported in Ref. [^. In the SM 
expectation of 0s (^, additional contributions to the leading b —)■ ccs tree Feynman 
diagram, as shown in Fig. are assumed to be negligible. However, the shift in 0s due to 
these contributions, called hereafter “penguin pollution”, is difficult to compute due to 
the non-perturbative nature of the quantum chromodynamics (QCD) processes involved. 
This penguin pollution must be measured or limited before using the 0s measurement 
in searches for BSM effects, since a shift in this phase caused by penguin diagrams is 
possible. Various methods to address this problem have been proposed j9-14 , and LHCb 
has recently published upper limits on the size of the penguin-induced phase shift using 
B^ —)■ J/0p° decays 


15 


Tree and penguin diagrams contributing to both B^ —>■ J/00 and B^ —>■ J/'ijjK'' 


*0 


decays are shown in Fig. | 
using —>■ J/'tpK*^ decay: 


first proposed in Ref. for the B^ 
Bg —)■ J/'ijjK*^ decay in Refs, [h 


In this paper, the penguin pollution in 0^ is investigated 
with J/0 —)■ and K*^ —)■ K~ti^, following the method 
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—)■ decay and later also discussed for the 

This approach requires the measurement of the 
branching fraction, direct CP asymmetries, and polarisation fractions of the —)■ J/ipK*^ 

decay. The measurements use data from proton-proton (pp) collisions recorded with the 
LHCb detector corresponding to 3.0 fb“^ of integrated luminosity, of which 1.0 (2.0) fb“^ was 
collected in 2011 (2012) at a centre-of-mass energy of 7 (8) TeV. The LHCb collaboration 
previously reported a measurement of the branching fraction and the polarisation fractions 
using data corresponding to 0.37fb“^ of integrated luminosity 16 . 


The paper is organised as follows: a description of the LHCb detector, reconstruction 


and simulation software is given in Sect. the selection of the B^ —)■ signal 

candidates and the —)■ J/ipK*^ control channel are presented in Sect. and the 

treatment of background in Sect. The invariant mass fit is detailed in Sect. 

The angular analysis and CP asymmetry measurements, both performed on weighted 
distributions where the background is statistically subtracted using the sPlot technique 
are detailed in Sect. The measurement of the branching fraction is explained in Sect. 
The evaluation of systematic uncertainties is described in Sect. along with the results, 
and in Sect, [^constraints on the penguin pollution are evaluated and discussed. 


Charge conjugation is implicit throughout this paper, unless otherwise specihed. 
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Figure 1: Decay topologies contributing to the Bg —>■ J/ipcj) channel (a,b) and —>■ 

channel (c,d). The tree diagrams (a,c) are shown on the left and the penguin diagrams (b,d) on 
the right. 


2 Experimental setup 

The LHCb detector [Isp^ is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector sur¬ 
rounding the pp interaction region, a large-area silicon-strip detector located upstream of 
a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The minimum distance 
of a track to a primary vertex, the impact parameter, is measured with a resolution of 
(15-|-29/pt) where pt is the component of the momentum transverse to the beam, 
in GeV/c. Different types of charged hadrons are distinguished using information from 
two ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a 
calorimeter system consisting of scintillating-pad and preshower detectors, an electromag- 
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netic calorimeter and a hadronic calorimeter. Muons are identified by a system composed 
of alternating layers of iron and multiwire proportional chambers. 

The online event selection is performed by a trigger, which consists of a hardware 
stage, based on information from the calorimeter and muon systems, followed by a 
software stage, which applies a full event reconstruction. In this analysis, candidates 
are first required to pass the hardware trigger, which selects muons with a transverse 
momentum px > 1.48GeV/c in the 7TeV data or px > 1.76GeV/c in the 8TeV data. In 
the subsequent software trigger, at least one of the final-state particles is required to have 
both px > 0.8 GeV/c and impact parameter larger than 100 pm with respect to all of the 
primary pp interaction vertices (PVs) in the event. Finally, the tracks of two or more of 
the final-state particles are required to form a vertex that is significantly displaced from 
any PV. Further selection requirements are applied offline in order to increase the signal 
purity. 

In the simulation, pp collisions are generated using Pythia 20,21 with a specific 


LHGb configuration |^. Decays of hadronic particles are described by EvtGen [23] , 
in which final-state radiation is generated using Photos 24 . The interaction of the 


generated particles with the detector, and its response, are implemented using the Geant 4 
toolkit 


25,26 as described in Ref. 27 


3 Event selection 


The selection of —?■ J/V’ K*^ candidates consists of two steps: a preselection consisting 

of discrete cuts, followed by a specific requirement on a boosted decision tree with gradient 
boosting (BDTG) [^,29 to suppress combinatorial background. All charged particles 
are required to have a transverse momentum in excess of 0.5 GeV/c^ and to be positively 
identified as muons, kaons or pions. The tracks are fitted to a common vertex which is 
required to be of good quality and significantly displaced from any PV in the event. The 
flight direction can be described as a vector between the production and decay vertices; 
the cosine of the angle between this vector and the 5° momentum vector is required to be 
greater than 0.999. Reconstructed invariant masses of the J/'ifj and K*^ candidates are 
required to be in the ranges 2947 < m^+^- < 3247 MeV/c^ and 826 < mK-n+ < 966 MeV/c^. 
The invariant mass is reconstructed by constraining the J/'ip candidate to its nominal 
, and is required to be in the range 5150 < < 5650MeV/c^. 


mass 
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The training of the BDTG is performed independently for 2011 and 2012 data, using 
information from the candidates: time of flight, transverse momentum, impact pa¬ 
rameter with respect to the production vertex and of fho decay vertex fit. The data 
sample used to train the BDTG uses less stringent particle identification requirements. 
When training the BDTG, simulated B^ —)■ J/'ipK*^ events are used to represent the signal, 
while candidates reconstructed from data events with J/'ipK~TT^ invariant mass above 
5401 MeV/c^ are used to represent the background. The optimal threshold for the BDTG 
is chosen independently for 2011 and 2012 data and maximises the effective signal yield. 
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Table 1: Expected yields of each background component in the signal mass range. 


Background sources 2011 data 2012 data 


B^ 

-)■ 

J/'07r+7r 

51 

± 

10 

115 

± 

23 


-)■ 

J/'07r+7r“ 

9.3 

± 

2.1 

25.0 

± 

5.4 

B^s 

-)■ 

J/^jK+K- 

10.1 

± 

2.3 

19.2 

± 

4.0 


-)■ 

J/i\)pK- 

36 

± 

17 

90 

± 

43 


-)■ 

JllppT^~ 

13.8 

± 

5.3 

27.3 

± 

9.0 


4 Treatment of peaking backgrounds 

After the suppression of most background with particle identification criteria, simulations 
show residual contributions from the backgronnds —)■ —)■ J/ip K~^K~, 

Bg —)■ J/ip'K^n~, and —)■ J/'07r’''7r“. The invariant mass distribntions of misidentihed 
B^ —)■ J/ip'K^n~ and B^ —)■ events peak near the 5° —)■ J/'ipK~'K^ signal peak 

dne to the effect of a wrong-mass hypothesis, and the misidentihed B^ —)■ J/'ijjK~^K~ 
candidates are located in the vicinity of the B^ —)■ J/'il)K^7r~ signal peak. It is therefore 
not possible to separate snch backgronnd from signal nsing information based solely on 
the invariant mass of the system. Moreover the shape of the rehected invariant 

mass distribution is sensitive to the daughter particles momenta. Due to these correlations 
it is difficult to add the 6-hadron to (where h is either a pion, a kaon or a proton) 

misidentihed backgrounds as extra modes to the ht to the invariant mass distribntion. 
Instead, simulated events are added to the data sample with negative weights in order to 
cancel the contribntion from those peaking backgronnds, as done previonsly in Ref. |^. 
Simulated 6-hadron to events are generated using a phase-space model, and then 

weighted on an event-by-event basis nsing the latest amplitnde analyses of the decays 

, and —)■ J/xljir^TT~ 

The snm of weights of each decay mode is normalised snch that the injected simnlated 
events cancel ont the expected yield in data of the specihc backgronnd decay mode. 

In addition to —)■ J/'tppK~ and B —)■ J/tl)h'^h~ decays, backgronnd from A^ —)• 
J/'ipp'n~ is also expected. However, in Ref. a full amplitude analysis was not performed. 
For this reason, as well as the fact that the A^ decays have broad mass distributions, the 
contribution is explicitly included in the mass ht described in the next section. Expected 
yields for both B —>■ J/'iph'^h~ and A^ —)■ J/'il)ph~ background decays are given in Table 
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Al J/^pK- pH, RO ^ J/ijK+K- pg, RO ^ J/V^Tr+Tr- 33; 


5 Fit to the invariant mass distribution 

After adding simulated R° —)■ J/^tt+tt", R° — J/^’ tt+tt” , R° — J/'ipK^K~, 
and A^ —)■ J/'il)pK~ events with negative weights, the remaining sample consists of 
R° —)■ J/ijj K~^7r~, Bg —>■ J/ijj K~7 i~^, A^ — >■ J/ippTT~ decays, and combinatorial background. 
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These four modes are statistically disentangled through a ht to the J/'ijj invariant 

mass. The combinatorial background is described by an exponential distribution, the 
—)■ J/'iljpTi~ decay by the Amoroso distribution 36 and the and 5° signals by the 
donble-sided Hypatia distribntion 1^, 


cr, A, C,/?, 01 , 02 , ni,n 2 ) oc 

A 
C 


((m — pY + e/3(m-n)^ 




li m — p < —Oi(j , 
if m — /r > 02(7, 
(^<y^y {m — pY + 6^^ otherwise. 


( 1 ) 


where Ky{z) is the modihed Bessel function of the second kind, 5 


a 


Kx+iiC)- 


a = ’’ A,B,C,D are obtained by imposing continuity and differentia¬ 

bility. This function is chosen because the event-by-event uncertainty on the mass has 
a dependence on the particle momenta. The estimate of the nnmber of B^ -A- J/tl)K^n~ 
decays lying under the peak is very sensitive to the modelling of the tails of the B^ 
peak. The fitted fraction is in good agreement with the estimate from simnlation. 

In the fit to data, the mean and resolution parameters of both the 5° and B^ Hypatia 
fnnctions are free to vary. All the remaining parameters, namely A, ai, ui, 02 and n 2 , are 
fixed to values determined from fits to 5° and B^ simulated events. All the A° —)■ J/'^pvr” 
shape parameters are fixed to values obtained from fits to simulated -A- J/V’p7r“ events, 
while the exponent of the combinatorial background is free to vary. 

Due to the small expected yield of A° —)■ J/Yp'n'~ decays compared to those of the 
other modes determined in the fit to data, and to the broad distribntion of A^ -A- J/'iljp'^~ 
decays across the J/'il>K~7r^ invariant mass spectrum, its yield is included in the fit as a 
Ganssian constraint nsing the expected nnmber of events and its nncertainties, as shown 
in Table [U 

From stndies of simnlated (MC) samples, it is fonnd that the resolntion of B^ and B^ 
mass peaks depends on both mK-n+ and cos(6*^), where 6^ is one of the helicity angles 
nsed in the angnlar analysis as defined in Sect. The fit to the J/YK~n^ invariant mass 
spectrum, including the evaluation of the sWeights, is performed separately in twenty 
bins, corresponding to four mK--K+ bins of 35MeV/c^ width, and hve equal bins in cos(6'^). 
The overall B^ and B^ yields are obtained from the sum of yields in the twenty bins, giving 


Nbo = 208656 ± 462 (stat)l^® (syst), (2) 

Nbo = 1808 ± 51 (stat)l 33 (syst), (3) 


where the statistical uncertainties are obtained from the quadratic sum of the uncertainties 
determined in each of the individual fits. Systematic uncertainties are discussed in Sect. 
The correlation between the and B^ yields in each bin are found to be smaller than 4%. 
The ratio of the and B^ yields is found to be Nbo/Nbo = (8.66±0.24(stat)^Q;}g (syst)) x 
10“^. Figure shows the sum of the ht results for each bin, overlaid on the 
mass spectrnm for the selected data sample. 
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Figure 2: The J/ijj invariant mass distribution with the sum of the fit projections in the 20 

mj^-^+ and cos{9^) bins. Points with error bars show the data. The projection of the fit result 
is represented by the solid blue line, and the contributions from the different components are 
detailed in the legend. At this scale the contribution of the —)• J/'^pvr” is barely visible. All 
the other peaking background components are subtracted as described in the text. 


6 Angular analysis 


6.1 Angular formalism 

This analysis uses the decay angles dehned in the helicity basis. The helicity angles are 
denoted by {9k, 9^, ^ph), as shown in Fig. The polar angle 9k {9^) is the angle between 
the kaon (/i+) momentum and the direction opposite to the momentum in the K~'k^ 
{pApr) centre-of-mass system. The azimuthal angle between the K~tt^ and pApT decay 
planes is (ph- The dehnitions are the same for or decays. They are also the same 
for —)■ decays. 

The shape of the angular distribution of —)■ J/'ipK*^ decays is given by Ref. 
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<Ar{9K,9f,,p>h) 


oc 


E 


|A|<J 

E 

A.J 


2J + 1 
dvr 




(4) 


where A = 0, ±1 is the J/'ip helicity, = ±1 is the helicity difference between the muons, 
J is the spin of the K~'k^ system, "H are the helicity amplitudes, and d are the small 
Wigner matrices. 

The helicity amplitudes are rotated into transversity amplitudes, which correspond to 
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final P eigenstates, 


.4s=WS, 

(5) 

^ = Wj, 

(6) 


(7) 


(8) 


The distribution in Eq. |^can be written as the sum of ten angular terms, four corresponding 
to the square of the transversity amplitude of each final state polarisation, and six 
corresponding to the cross terms describing interference among the final polarisations. 

The modulus of a given transversity amplitude, A^, is written as \Ax\, and its phase as 
5x- The convention 5o = 0 is used in this paper. The P-wave polarisation fractions are 
fi = l^ip/dydop + l^iiP + with i = 0, II, T and the S-wave fraction is defined as 

Fs = l^sp/dAoP + |A||p + |v4_i_d + |v4gp). The distribution of the CP-conjugate decay 
is obtained by flipping the sign of the interference terms which contain |A_l|. For the 
CP-conjugate case, the amplitudes are denoted as Ai. Each Ai and the corresponding Ai 


are related through the CP asymmetries, as described in Sect. 6.3 


6.2 Partial-wave interference factors 

In the general case, the transversity amplitudes of the angular model depend on the K~'k^ 
mass {mK--K+)- This variable is limited to be inside a window of ±70MeV/C around the 
K*^ mass. Figure]^ shows the efficiency-corrected mK-n+ spectra for P° and P° using the 
nominal sets of sWeights. 

In order to account for the mK-n+ dependence while keeping the framework of an 
angular-only analysis, a fit is performed simultaneously in the same four mK--K+ bins 
defined in Sect. Different values of the parameters |^sP and are allowed for each 
bin, but the angular distribution still contains mass-dependent terms associated with the 
interference between partial-waves. If only the S-wave and P-wave are considered, such 
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interference terms correspond to the following complex integrals, 


V X S* ^ EmirriKn) drriK-K 


\P\^^ emirriKn) drriKn JZl" $ emimK-n) drriKn 


H 

Kt\ 

Ktt 


— Cspe 


-idsp 


(9) 


where lower (higher) limit of the bin, SmifEiK-K) is the acceptance for a K~tt^ 

candidate with mass mx-K (see Appendixfor a discussion on the angular acceptance), 
<I> stands for the phase space, and V (5) is the P-wave (S-wave) propagator. The phase 
space term is computed as 


$ 


pq 


m 


2 ’ 
Ktt 


( 10 ) 


where p denotes the K*^ momentum in the rest frame and q refers to the K~ momentum 
in the K*^ rest frame. 

The phase ^sp is included in the dehnition of 6s but the Csp factors, corresponding 
to real numbers in the interval [0,1], have to be computed and input to the angular ht. 
The contribution of D-wave (J = 2) in the mx-TT+ range considered is expected to be 
negligible. Therefore the nominal model only includes S-wave and P-wave. To determine 
the systematic uncertainty due to possible D-wave contributions, Csd and Cpp factors 
are also computed, using analogous expressions to that given in Eq. The Cij factors are 
calculated by evaluating numerically the integrals using the propagators outlined below, 
and are included as fixed parameters in the ht. A systematic uncertainty associated to the 
different possible choices of the propagator models is afterwards evaluated. 


U 

> 




cd 

75 

a 

cd 

U 



"V-jr [MeV/c^] 


Figure 4: Efficiency corrected mx-Tr+ distribution for shown in squares (red) and shown 
in circles (black) using sWeights computed from the maximum likelihood fit to the 
invariant mass spectrum. 


The S-wave propagator is constructed using the LASS parametrisation , consisting 
of a linear combination of the Ag(1430)° resonance with a non-resonant term, coming 
from elastic scattering. The P-wave is described by a combination of the A*(892)° and 
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Table 2: The Csp, Csd and Cpd factors calculated in each of the four bins around the 

K*^ peak. 


Bin 

^k-tt+ range (MeV/c^) 

C'sp 

C'sD 

C'PD 

0 

[826, 861] 

0.968 ± 0.017 

0.9968 ± 0.0030 

0.9827 ± 0.0048 

1 

[861, 896] 

0.931 ± 0.012 

0.9978 ± 0.0021 

0.9402 ± 0.0048 

2 

[896, 931] 

0.952 ± 0.012 

0.9983 ± 0.0016 

0.9421 ± 0.0056 

3 

[931, 966] 

0.988 ± 0.011 

0.9986 ± 0.0012 

0.9802 ± 0.0066 


resonances using the isobar model 40 , and the D-wave is assumed to come from 
the 4 ^ 2 ( 1430 )° contribution. Relativistic Breit-Wigner functions, multiplied by angular 
momentum barrier factors, are used to parametrise the different resonances. Table 
contains the computed Csp, Csd and C-pD factors. 


6.3 CP asymmetries 

The direct CP violation asymmetry in the decay rate to the final state /(s)*, with 
fs,i = J/il^ {K~TT~^)i and fi = J/^lJ {K~^7i~)i, is dehned as 




0 V f '1 _ l^(^) 

(s) Jis)i) - yy 


|^(s) ) 


+ |^(s)) 


( 11 ) 


where R(s)i are the transversity amplitudes dehned in Sect. 6.1 and the additional index s 


is used to distinguish the and the R°-meson. The index i refers to the polarisation of 
the hnal state {i = 0, ||, T, S) and is dropped in the rest of this section, for clarity. 

The raw CP asymmetry is expressed in terms of the number of observed candidates by 


fis)) — 




( 12 ) 


Both asymmetries in Eq. 11 and Eq. are related by 

^ /(,,) ~ AZiB'l.^ ^ ft.)) - Ct.tMf) - 't(,)Ap(B»)), ( 13 ) 


where Rd(/) is the detection asymmetry, dehned as in Eq, 
production asymmetry, dehned as in Ec 


. 3lp(B».,) is the B»,)-By 


dilution due to 5°^^—5°^^ oscillations 


42' 


, (15), C(s) = +1(—1) and K(s) accounts for the 
. The K(s) factor is evaluated by 


K-{s) — 




r(^)tcosh(^^ 


fo e 


t]s(t)dt 


(14) 


where e(t) is the time-dependent acceptance function, assumed to be identical for the 
Bg —)■ J/iIjK*^ and B^ —)■ J/iIjK*^ decays. The symbols r( 5 ) and Am(s) denote the decay 
width and mass differences between the mass eigenstates. 
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The production asymmetry is defined as 


(S?,|) 


(g?.)) - (fl(.)) 

( 5 ?.)) + («(.)) ’ 


(15) 


where a is the production cross-section within the LHCb acceptance. The production 


asymmetries reported in Ref. 
obtain 
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are reweighted in bins of transverse momentum to 


Ap{B^) = (—1.04 ± 0.48 (stat) ± 0.14 (syst)) % , 
ylp(R°) = (-1.64 ± 2.28 (stat) ± 0.55 (syst)) % . 


The K(s) factor in Eq. 14 is determined by fixing Ar(s), Am(s) and r(s) to their world 
average values and by fitting the decay time acceptance e{t) to the nominal data 
sample after applying the B^ sWeights, in a similar way to Ref. |^. It is equal to 0.06% 
for B^ decays, and 41% for B^. This reduces the effect of the production asymmetries to 
the level of 10“® for —)■ and 10“^ for B^ —)■ decays. 

Other sources of asymmetries arise from the different final-state particle interactions 
with the detector, event reconstruction and detector acceptance. The detection asymmetry, 
A-D^f), is defined in terms of the detection efficiencies of the final states, 


as 


Mf) 


gdet(y) _ 
^det(J)+^det(J) ■ 


(16) 


The detection asymmetry, measured in bins of the momentum in Ref. [^, is weighted 
with the momentum distribution of the kaon from the —?• J/ip decays to give 


At){B^) = ( 1.12 ±0.55 (stat)) %, 
ylD(5°) = (-1.09 ±0.53 (stat)) %. 


7 Measurement of 

The branching fraction B{Bg —)• J/'ipK*^) is obtained by normalising to two different 
channels, B^ —)■ J/V’0 and B^ —)■ J/'ipK*^, and then averaging the results. The expression 

g(R0 ^ X ^ K+n-) ^ N^q^j/^k^q U 

B{Bq —)■ J/i>X) X B{X — )■ h+h-) Xsq^j/ipx ^BO^j/i>K*o fs ’ 

is used for the normalisation to a given Bq —)■ J/^l)X decay, where N refers to the yield of 
the given decay, £ corresponds to the total (reconstruction, trigger and selection) efficiency, 
and fq = fsifd) are the R°(i?°)-meson hadronisation fractions. 
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The event selection of —)■ J/ipcj) candidates consists of the same reqnirements as 

those for —)■ J/ipK*^ candidates (see Sect. [^, with the exception that 0 candidates are 
reconstructed in the K^K~ state so there are no pions among the final state particles. In 
addition to the other requirements, reconstructed 0 candidates are required to have mass 
in the range 1000 < mK-K+ < 1040MeV/c^ and to have a transverse momentum in excess 
of IGeV/cT 


7.1 Efficiencies obtained in simulation 


A first estimate of the efficiency ratios is taken from simulated events, where the particle 


identihcation variables are calibrated using D*^ decays. The efficiency ratios estimated from 
simulation, for 2011 (2012) data, are = 0.929T0.012 (0.927T0.012) 

and = 1-991 ± 0.025 (1.986 ±0.027). 


7.2 Correction factors for yields and efficiencies 


The signal and normalisation channel yields obtained from a mass fit are affected by the 
presence of a non-resonant S-wave background as well as interference between S-wave and 
P-wave components. Such interference would integrate to zero for a flat angular acceptance, 
but not for experimental data that are subject to an angle-dependent acceptance. In 
addition, the efficiencies determined in simulation correspond to events generated with an 
angular distribution different from that in data; therefore the angular integrated efficiency 
also needs to be modified with respect to simulation estimates. These effects are taken 
into account using a correction factor cu, which is the product of the correction factor to 
the angular-integrated efficiency and the correction factor to the P-wave yield: 


AT — AT — c^-lVlL; 

^ ^ TvlC ^ ^ 

NBq^J/ipX ^B°^J/i)K*° NBq^J/tpx ^B0^J/ipK*0 ^B^^J/'ipK*^ 


(18) 


where Nsq^j/^px are the yields obtained from the mass fits, 

are the efficiencies obtained in simulation, and u is calculated as 


^Bq^J/lPX 


frx 

^ Bq^JjilX 
CBq-^JIi)X 


(19) 


where is the fraction of the P-wave X resonance in a given Bq —)• 7/0 X decay 

(related to the presence of S-wave and its interference with the P-wave), and CSq^j/iijx is 
a correction to simulated values of the decay parameters 

differ slightly from those measured. The values obtained for the uj correction factors are 


^bo^jIiIjK*o = 1-149 ± 0.044 (stat) ± 0.018 (syst), 
UBO^j/^jjK*o = 1-107 ± 0.003 (stat) ± 0.038 (syst), 
uiBo^j/tp<l> = 1-013 ± 0.002 (stat) ± 0.007 (syst). 
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7.3 Normalisation to 

The study of penguin pollution requires the calculation of ratios of absolute ampli¬ 
tudes between 5° —)■ and 5° —)■ J/V’0. Thus, normalising —)• J/ipK*^) to 

B{B^g —>■ J/ipcp) is very useful. This normalisation is given by 


^ Nbo^j/^k--k+ ^ ^ ^ B{(t>^K+K ) 

—)> J/b’0) ^BO^J/ipK+K- j/ipK*° B(K*^ —> K 7r+) 
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Using 


where B{K*^ = 2/3 and B{^ K+R-) = (49.5 ± 0.5)% 

NB0^j/-tjjK--K+ as given in Eq. and Nbo^j/^k+k- = 58 09f ± 243 (stat) ± 3f9 (syst) 
as obtained from a ht to the invariant mass of selected J/^0 candidates, where 

the signal is described by a double-sided Hypatia distribution and the combinatorial 
background is described by an exponential distribution, a value of 


B{B« ^ J/V^) = b05 ± 0.19(stat) ± 0.13(syst))% 

is obtained. 


7.4 Normalisation to B^ — 

The normalisation to B^ —)■ J/ipR*^ is given by 


J/tpR*^) _ NB0^j/^jiK--K+ fd ^B0^J/',pK*0 ^BO^Jf-ipK^o 

8{B» ^ Jf4,K>«) ^ '"s'' ’ 


where A^B 0 ej/yA-+ 7 r- and Nb^^j/^k-it*- are given in Eq. 0 and Eq. respectively, and 


^BO^J/^K*o 

^BO^J/ipK*o 


0.963 ± 0.036 (stat) ± 0.031 (syst), 


resulting in a value of 


ro v 

B{B^ ^ ^ ^ ^ ^ 

where the third uncertainty comes from the hadronisation fraction ratio 
/rf//, = 3.86 ± 0.22 0. 


7.5 Computation of B{B^ — 

By multiplying the fraction given in Eq. by the branching fraction of the decay 
B^ —)■ J/'ipR*^ measured at Bell^ (1.29 ± 0.05 (stat) ± 0.13 (syst)) x 10“^ 46 , and taking 


^The result from Belle was chosen rather than the PDG average, since it is the only B{B^ 
measurement that subtracts S-wave contributions. 
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into account the difference in production rates for the B^B and B^B^ pairs at the T{4S) 
resonance, i.e. T{B~^B~)/T{B^B^) = 1.058 ± 0.024 0 , the value 

^ J/^K*^)d = (3.95 ±0.18 (stat) ± 0.16 (syst) ± 0.23 {fd/fs) 

±0.43 {B{B^ J/ipK*^))) X 10-® 

is obtained, where the fourth uncertainty arises from B{B^ —)■ J/'ipK*^). A second estimate 
of this quantity is found via the normalisation to —)■ J/'ip<p) |^, updated with the 

value of fd/fs from Ref. to give B{B'/ —)■ J/fjcj)) = (1.038 ± 0.013 (stat) ± 0.063 (syst) ± 
0.060 {fd/fs)) X 10“^, resulting in a value of 

B{Bl J/i/K*^)^ = (4.20 ± 0.20 (stat) ± 0.13 (syst) ± 0.36 (i3(R° ^ ^/V'^))) x lO”*^, 

where the third uncertainty comes from B{B^^ —)■ J/ipf). Both values are compatible within 
uncorrelated systematic uncertainties and are combined, taking account of correlations, to 
give 


i5(R° ^ J/fjK*^) = (4.14 ± 0.18 (stat) ± 0.26 (syst) ± 0.24 (fd/fs)) x lO”® , 
which is in good agreement with the previous LHCb measurement 16 , of (4.41 q 4 ± 0.8) x 
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8 Results and systematic uncertainties 


Section |8.1| presents the results of the angular ht as well as the procedure used to estimate 
the systematic uncertainties, while in Sect. | 8 . 2 | the results of the branching fraction 


measurements and the corresponding estimated systematic uncertainties are discussed. 


8.1 Angular parameters and CP asymmetries 

The results obtained from the angular ht to the B^ —)■ J/tf) K*^ events are given in Table 
and Table for the P-wave and S-wave parameters, respectively. For comparison, the 
previous LHCb measurements of /o and f\\ were 0.50 ± 0.08 ± 0.02 and 0.19 ](q;o 8 ± 0.02, 
respectively. The angular distribution of the signal and the projection of the htted 
distribution are shown in Fig. The statistical-only correlation matrix as obtained from 
the ht to data is given in Appendix The polarisation-dependent CP asymmetries are 
compatible with zero, as expected in the SM. The polarisation fractions are in good 


agreement with the previous measurements 16 performed on the same decay mode by the 
LHCb collaboration using data corresponding to an integrated luminosity of 0.37fb~^. 

Various sources of systematic uncertainties on the parameters of the angular ht are 
studied, as summarised in Table and Table for the P-wave and S-wave parameters. 
Two classes of systematic uncertainties are dehned, one from the angular ht model and 
another from the mass ht model. Since the angular ht is performed on the data weighted 
using the signal sWeights calculated from the ht to the J/fjK~n^ invariant mass, biases on 
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the mass fit results may be propagated to the sWeights and thus to the angular parameters. 
Overall, two sources of systematic uncertainties dominate: the angular acceptance and the 
correlation between the J/'ipK~n~^ invariant mass and 9^. 

8.1.1 Systematic uncertainties related to the mass fit model 

To determine the systematic uncertainty arising from the hxed parameters in the description 
of the J/tj) K~tt^ invariant mass, these parameters are varied inside their uncertainties, as 
determined from hts to simulated events. The ht is then repeated and the widths of the 
and yield distributions are taken as systematic uncertainties on the branching fractions. 
Correlations among the parameters obtained from simulation are taken into account in 
this procedure. For each new ht to the J/'il)K~7r~^ invariant mass, the corresponding set of 
sWeights is calculated and the ht to the weighted angular distributions is repeated. The 
widths of the distributions are taken as systematic uncertainties on the angular parameters. 
In addition, a systematic uncertainty is added to account for imperfections in the modelling 





(p^ (rad) 


-- S-P wave interference 

- S-wave 

. P-wave odd + P-P wave interference 

- P-wave even 

- Totai 


Figure 5: Fitted signal distributions compared with the weighted angular distributions with 
sWeights. Points with error bars show the data. The projection of the fit result is represented by 
the solid black line, and the contributions from the different amplitude components are described 
in the legend. 
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Table 3: Summary of the measured J/ipK*^ P-wave properties and their statistical and 

systematic uncertainties. When no value is given, it means an uncertainty below 5 x 10“^, except 
for the two phases, (fii (rad) and (f_L (rad), in which case the uncertainty is below 5 x 10“^. 



/o 

h 

5|l 



Af 

A<IP 

Fitted value 

0.497 

0.179 

-2.70 

0.01 

-0.048 

0.171 

-0.049 

Statistical uncertainties 

0.025 

0.027 

0.16 

0.11 

0.057 

0.152 

0.096 

Angular acceptance 

(simulation statistics) 

0.018 

0.008 

0.02 

0.01 

0.009 

0.017 

0.008 

Angular acceptance 

(data-simulation differences) 

0.015 

0.007 

0.17 

0.10 

0.007 

— 

0.015 

(PsP factors 

— 

0.001 

— 

— 

0.001 

0.002 

0.002 

D-wave contribution 

0.004 

0.003 

— 

— 

0.002 

0.015 

0.002 

Background 
angular model 

0.004 

0.002 

0.02 

0.01 

0.004 

+0.012 

-0.004 

0.002 

Mass parameters and 
contamination 

— 

— 

— 

— 

0.001 

0.001 

— 

Mass-cos(0^) 

correlations 

0.007 

0.006 

0.07 

+0.02 

-0.04 

0.014 

+0.009 

-0.012 

0.016 

Fit bias 

— 

0.001 

0.01 

0.07 

0.003 

0.002 

0.005 

Detection 

asymmetry 

Production 

— 

— 

— 

— 

0.005 

0.005 

0.006 

asymmetry 








Quadratic sum of 
systematic uncertainties 

0.025 

0.013 

0.19 

+0.012 

-0.013 

0.020 

+0.028 

-0.027 

0.025 

Total uncertainties 

0.035 

0.030 

0.25 

+0.016 

-0.017 

0.060 

0.154 

0.099 


of the upper tail of the and peaks. Indeed, in the Hypatia distribution model, the 
parameters 02 and n 2 take into account effects such as decays in flight of the hadron, that 
affect the lineshape of the upper tail and could modify the B^ leakage into the B^ peak. 
The estimate of this leakage is recalculated for extreme values of those parameters, and 
the maximum spread is conservatively added as a systematic uncertainty. 

Systematic uncertainties due to the hxed yields of the B^ —)■ J/'ip K~^K~, 
Bg —)■ J/ipTT^TT~, B^ —)■ J/07r+7r“, and —)■ J/%l:pK~ peaking backgrounds]^ are evalu- 

^The yields of the subtracted backgrounds can be considered as fixed, since the sum of negative weights 
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Table 4: Summary of the measured —)• J/ipK*^ S-wave properties and their statistical and 
systematic uncertainties. When no value is given, it means an uncertainty below 5 x 10“^, except 
for the four phases related to the S-wave component, (5s (I'ad), in which case the uncertainty is 
below 5 X 10“^. The binning definition is identical to the one given in Table ^ 




,binO 


„binl 

_bin2 





IIIr- 

'77+ 

IILr- 

'77+ 

IILr- 

'77 + 



Ts 


Fs 

<5s 

Fs 


Fs 

<5s 

Fitted value 

0.167 

0.475 

0.54 

0.080 

-0.53 

0.044 

-1.46 

0.523 

-1.76 

Statistical uncertainties 

0.114 

+0.108 

-0.112 

0.16 

+0.031 

-0.025 

+0.25 

-0.21 

+0.042 

-0.029 

+0.22 

-0.19 

+0.109 

-0.112 

+0.13 

-0.14 


Angular acceptance 

(simulation statistics) 

0.028 

0.039 

0.03 

0.012 

0.065 

0.015 

0.10 

0.065 

0.06 

Angular acceptance 

(data—simulation differences) 

0.015 

0.058 

0.08 

0.019 

0.18 

0.027 

0.27 

0.006 

0.04 

Csp factors 


0.002 

0.01 

0.001 

— 

0.002 

— 

0.001 

0.01 

D-wave contribution 

0.008 

0.010 

0.02 

0.005 

0.03 

0.008 

0.08 

0.002 

0.04 

Background 
angular model 

0.001 

0.002 

0.01 

+0.000 

-0.001 

0.01 

— 

0.03 

+0.002 

-0.000 

+0.07 

-0.04 

Mass parameters and 
contamination 

0.001 

0.001 

0.01 

— 

— 

— 

— 

— 

— 

Mass-cos (+) 
correlations 

+0.023 

-0.029 

+0.040 

-0.028 

0.05 

0.003 

0.04 

+0.006 

-0.016 

0.02 

+0.009 

-0.011 

0.03 

Fit bias 

0.004 

0.005 

0.01 

0.003 

0.02 

0.007 

0.032 

0.015 

0.01 

Detection 

asymmetry 

Production 

0.005 

— 

— 

— 

— 

— 

— 

— 

— 


asymmetry 











Quadratic sum of 
systematic uncertainties 

+0.041 

-0.044 

+0.081 

-0.076 

0.10 

0.023 

0.20 

+0.033 

-0.036 

0.30 

0.068 

+0.11 

-0.09 

Total uncertainties 

+0.120 

-0.122 

0.135 

0.19 

+0.039 

-0.034 

+0.32 

-0.29 

+0.054 

-0.047 

+0.37 

-0.35 

+0.128 

-0.131 

0.17 


ated by repeating the fit to the invariant mass varying the normalisation of all background 
sources by either plus or minus one standard deviation of its estimated yield. For each of 
the new mass hts, the angular ht is repeated using the corresponding new sets of sWeights. 
The deviations on each of the angular parameters are then added in quadrature. 

Correlations between the invariant mass and the cosine of the helicity angle 

9^ are taken into account in the nominal ht model, where the mass ht is performed in hve 
bins of cos(6*^). In order to evaluate systematic uncertainties due to these correlations, 
the mass ht is repeated with the full range of cos(6'^) divided into four or six equal bins. 

used to subtract them is constant in the nominal fit. 
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For each new mass fit, the angular ht is repeated using the corresponding set of sWeights. 
The deviations from the nominal result for each of the variations are summed quadratically 
and taken as the systematic uncertainty. 


8.1.2 Systematic uncertainties related to the angular fit model 


In order to account for systematic uncertainties due to the angular acceptance, two distinct 
effects are considered, as in Ref. [^. The first is due to the limited size of the simulation 
sample used in the acceptance estimation. It is estimated by varying the normalisation 
weights 200 times following a Gaussian distribution within a hve standard deviation range 
taking into account their correlations. For each of these sets of normalisation weights, 
the angular ht is repeated, resulting in a distribution for each htted parameter. The 
width of the resulting parameter distribution is taken as the systematic uncertainty. Note 
that in this procedure, the normalisation weights are varied independently in each mK-n+ 
bin. The second effect, labelled as data-simulation corrections in the tables, accounts for 
differences between the data and the simulation, using normalisation weights that are 
determined assuming the amplitudes measured in Ref. 47 . The difference with respect to 


the nominal ht is assigned as a systematic uncertainty. The uncertainties due to the choice 
of model for the Gsp factors are evaluated as the maximum diherences observed in the 
measured parameters when computing the Gsp factors with all of the alternative models, 
as discussed below. Instead of the nominal propagator for the S-wave, a combination of 
the iFQ(800)° and iFQ(1430)‘^ resonances with a non-resonant term using the isobar model 
is considered, as well as a K-matrix 48 version. A pure phase space term is also used, in 


order to account for the simplest possible parametrisation. For the P-wave, the alternative 
propagators considered are the K*{892)^ alone and a combination of this contribution 
with the iFj"(1410)° and the iFj"(1430)° using the isobar model. 

In order to account for the absence of D-wave terms in the nominal ht model a new ht 
is performed, including a D-wave component, where the related parameters are hxed to 
the values measured in the region. The differences in the measured parameters 

between the results obtained with and without a D-wave component are taken as the 
corresponding systematic uncertainty. 

The presence of biases in the ht model itself is studied using parametric simulation. For 
this study, 1000 pseudoexperiments were generated and htted using the nominal shapes, 
where the generated parameter values correspond to the ones obtained in the ht to data. 
The difference between the generated value and the mean of the distribution of htted 
parameter values are treated as a source of systematic uncertainty. 

Finally, the systematic uncertainties due to the hxed values of the detection and 
production asymmetries are estimated by varying their values by ±1 standard deviation 
and repeating the ht. 
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8.2 Branching fraction 


Several sources of systematic uncertainties on the branching fraction measurements are 
studied, summarised along with the results in Table systematic uncertainties due to the 
external parameter fd/fs and due to the branching fraction B{(j) —)■ K~^K~); systematic 
uncertainties due to the ratio of efficiencies obtained from simulation and due to the angular 
parameters, propagated into the oj factors (see Sect. 8.1); and systematic uncertainties 
affecting the J/tjj K*^ and —)■ J/tjj K*^ yields, which are determined from the ht to 

the J/'il>K^7r~ invariant mass and described in Sect. 8 
due to the B^ —)■ J/ijjcj) yield determined from the ht to the J/ijjK^K 


Finally, a systematic uncertainty 
invariant mass 


distribution, described in Sect. |7.3[ is also taken into account, where only the effect due 
to the modelling of the upper tail of the B^ peak is considered (see Sect, 
computation of the absolute branching fraction i3(i?° —)■ J/'ipK*^) (see 
additional systematic sources are taken into account, the uncertainties in the external 
parameters i3(il° —>■ J/'ijjK*^) and B{B^g — J/ixf)). 


00 

. For the 

Sect. 

7.5 

), two 


Table 5: Summary of the measured values for the relative branching fractions and their statistical 
and systematic uncertainties. 


Relative branching fraction 


Nominal value 

2.99 

4.05 

Statistical uncertainties 

0.14 

0.19 

Efficiency ratio 

0.04 

0.05 

Angular correction (u) 

0.09 

0.07 

Mass model (effect on the yield) 

0.06 

0.08 

fd/fs 

0.17 

— 

B{(j) K+K-) 

— 

0.04 

Quadratic sum (excluding fd/fs) 

0.12 

0.13 

Total uncertainties 

0.25 

0.23 


9 Penguin pollution in (j)s 

9.1 Information from 

Following the strategy proposed in Refs. [^[^[^, the measured branching fraction, 
polarisation fractions and CP asymmetries can be used to quantify the contributions 
originating from the penguin topologies in —)■ To that end, the transition 

amplitude for the B^ —)■ J/'ipK*^ decay is written in the general form 

R (R° ^ = -AA [1 - , (23) 
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where A = \ Vus\ = 0.22548tQ;oQ034 and i labels the different polarisation states. In the 
above expression, Ai is a CP-conserving hadronic matrix element that represents the tree 
topology, and a* parametrises the relative contribntion from the pengnin topologies. The 
CP-conserving phase difference between the two terms is parametrised by 0*, whereas their 
weak phase difference is given by the angle 7 of the Unitarity Triangle. 

Both the branching fraction and the CP asymmetries depend on the penguin parameters 
Oj and Oi- The dependence of is given by 


Air = 


2 aj sin 6i sin 7 


2 aj cos 9i cos 7 + 


(24) 


To use the branching fraction information an observable is constructed |^: 


H, = - 


A' 


Ai 


$ 






$ 


^ —>■ J/ljj (j))theo f i 


7 ’ 

i 


^bO ’ ^bO ^ 

1 — 2 aj cos 9i cos 7 + 

1 -I- 2ea' cos 9'- cos 7 -|- e^af 

where fl''^ represents the polarisation fraction, 

A2 


e = 


1-A2 


= 0.0536 ± 0.0003 6 , 


(25) 


(26) 


and ^{x,y) = a/(1 - {x — |/)^)(1 — (a; -|- yy) is the standard two-body phase-space func¬ 
tion. The primed quantities refer to the P° —)■ J/ijj (j) channel, while the non-primed ones 
refer to P° —?■ J/'ipK*^. The penguin parameters a' and 9[ are dehned in analogy to a* and 
and parametrise the transition amplitude of the P° —?• decay as 


Al(P°^(J/^0),) 



1 -|- ea'e 


i 6 »' 



(27) 


Assuming SU(3) flavour symmetry, and neglecting contributions from exchange and 
penguin-annihilation topologies,]^ which are present in P° —>■ J /'00 but have no counterpart 
in P° — J/ipK*^, we can identify 


a, = a,- 


9'= 9,. 


(28) 


The contributions from the additional decay topologies in P° —)■ J/ilxj) can be probed 
using the decay P° —>■ (^. The current upper limit on its branching fraction is 

P(P° —)> J/V’0) < 1.9 X 10“^ at 90% conhdence level (C.L.) j^, which implies that the size 
of these additional contributions is small compared to those associated with the penguin 
topologies. 


"‘We follow the decomposition introduced in Ref. 


49 
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The Hi observables are constructed in terms of the theoretical branching fractions 
dehned at zero decay time, which differ from the measured time-integrated branching 
fractions 51 due to the non-zero decay-width difference AF^ of the meson system [^. 


The conversion factor between the two branching fraction definitions 51 is taken to be 

B{B —)■ /)theo 


Vs 


B{B f) 1 - ysTji cos(0f’^) 


(29) 


where rji is the CP eigenvalue of the final state, and i/s = Ars/2r^. Taking values for F^, 
AF^ and 0®’^ from Refs. (^[^, the conversion factor is 1.0608±0.0045 (0.9392±0.0045) for 
the CP-even (-odd) states. For the flavour-specific P° —)• J/'ipK*^ decay rji = 0, resulting 
in a conversion factor of 0.9963 ± 0.0006. The ratios of hadronic amplitudes are 


calculated in Ref. 52 following the method described in Ref. 53 


results on form factors from Light Cone QCD Sum Rules (LCSR) 


and using the latest 
This leads to 
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Ho = 0.98 ± 0.07 (stat) ± 0.06 (syst) ± 0.26 {\AJA\) , 

P|l = 0.90 ± 0.14 (stat) ± 0.08 (syst) ± 0.21 {\MJA\) , 

H^ = 1.46 ± 0.14 (stat) ± 0.11 (syst) ± 0.28 {\AJA\) ■ 


Assuming Eq. 


28 


and external input on the Unitarity Triangle angle 7 = ( 73 . 2 ]'^ 7 q)° 6 , 
the penguin parameters a, and 9i are obtained from a modihed least-squares ht to H^ 

in Eq. [^and Eq. 25 The information on 7 is included as a Gaussian constraint in the ht. 


The values obtained for the penguin parameters are 


Oo 

— n 04+0-95 

— U.UU_g Q4 , 

9o = 


a\\ 

— U.OZ_g 32 , 

9\\ = 

-(I5l}f)° 

a± 

— 0 44+0-^^ 

— U.UU_g 27 , 

9± = 



For the longitudinal polarisation state the phase 9 is unconstrained. Correlations between 
the experimental inputs are ignored, but the effect of including them is small. The 
two-dimensional conhdence level contours are given in Fig. This hgure also shows, 
as different (coloured) bands, the constraints on the penguin parameters derived from 
the individual observables entering the ht. The thick inner darker line represents the 
contour associated with the central value of the input quantity, while the outer darker 
lines represent the contours associated with the one standard deviation changes. For 
the parallel polarisation the central value of the H observable does not lead to physical 
solutions in the plane, and the thick inner line is thus absent. 

When decomposed into its different sources, the angle (ps takes the form 


-'S,! 




(30) 


where —2/3s is the SM contribution, 0^®^ is a possible BSM phase, and is a shift 

introduced by the presence of penguin pollution in the decay B^ —)• J/tpxp. In terms of the 
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Figure 6: Limits on the penguin parameters a* and 6i obtained from intersecting contours derived 
from the CP asymmetries and branching fraction information in J/ipK*^. Superimposed 

are the confidence level contours obtained from a fit to the data. Shown are the longitudinal 
(top), parallel (middle) and perpendicular (bottom) polarisation. 
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penguin parameters a' and 6 *', the shift is defined as 

.. , J/P0X ^ 2e< cos Q\ sin 7 + sin(27) 

1 + 2 ea' cos 6 *' cos 7 + e^af cos( 27 ) 


(31) 


Using Eqs. 
giving 


and 31, the ht results on a* and 6 *j given above constrain this phase shift, 


J/plf) _ 

0.003 

(stat) 

^aooQ (syst) 

to 0 

CO 

0 0 
d d 
+ 1 

Jlplj) _ 

«.ll 

0.031 

^^37 (stat) 

1S:S1“ (syst)±0.032 (|yl'/AI). 

J/plf) _ 

-0.045±0.012 (stat)±0.008 (syst) 



which is in good agreement with the values measured in Ref. (Tsl, and with the predictions 


given in Refs. 12 -14 


The above results are obtained assuming SU(3) flavour symmetry and neglecting con¬ 
tributions from additional decay topologies. Because Oje'®* represents a ratio of hadronic 
amplitudes, the leading factorisable SU(3)-breaking effects cancel, and the relation be¬ 
tween and is only affected by non-factorisable SU(3)-breaking. This can be 
parametrised using two SU(3)-breaking parameters ^ and 5 as 


a[ = ^ X tti, 9 [ = 9 i + 6 . 


(32) 


The above quoted results assume ^ = 1 and 5 = 0. The dependence of the uncertainty 
on on the uncertainty on ^ is illustrated in Fig. 7, while the dependence on the 

uncertainty on 6 is negligible for the solutions obtained to {ai,9i}. 


9.2 Combination with 

The information on the penguin parameters obtained from —>■ Jf-ipK*^ can be comple¬ 
mented with similar information from the SU(3)-related mode —>■ J/'iljp^ [^. Both 
modes describe a 6 —)■ ccd transition, and are related by exchanging the spectator s ^ d 
quarks. The decay amplitude of —)■ J/0 is also parametrised as 


A{B^ ^{J/ijp%) =-XA 




(33) 


which is the equivalent of Eq. 23 In contrast to B^ —)■ J/ipK*^, however, di and 9i also 


include contributions from exchange and penguin-annihilation topologies, which are present 
in B^ —)• J/'ipp^ but have no counterpart in —)■ Assuming SU(3) symmetry, 

and neglecting the contributions from the additional decay topologies in —)■ J/il^cj) and 
B^ —)■ J/ipp^, the relation in Eq. 28 can be extended to 


Qj^ - CL-j - Ctr} 


9'i = 9i = 9i 


(34) 


which allows a combined fit to be performed to the CP asymmetries and branching fraction 
information in B^ —)■ J/ipK*^ and B^ —)■ J/'ifjp^. 
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Figure 7: Dependence of the uncertainty on the penguin shift on the uncertainty on 

The bands correspond to the 68% C.L. The longitudinal (top), parallel (middle) and perpendicular 
(bottom) polarisations are shown. 
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The —)■ J/'ifj decay exhibits decay-time-dependeiit CP violation, which is described 

by two parameters, the direct CP asymmetry Cj, which in the SU(3) limit is related to 
as Ci = —Af^, and the mixing-induced CP asymmetry Si. Their dependence on the 
penguin parameters a* and 9i is given by 




-Vi 


2 tti sin 6 i sin 7 
1 — 2 cij cos 6 i cos 7 + ’ 

sin (j)d- 2 di cos 6 i sm{(j)d + 7) -f- sin(0d + 27) 
1 — 2 Oj cos 9i cos 7 + 


(35) 

(36) 


where rji is the polarisation-dependent CP eigenvalue of the —)■ J/tp decay, and pd 

is a CP-violating phase arising from the interference between B^-B^ mixing and the 
subsequent B^ decay. The use of Si to constrain the penguin parameters a* and 9i requires 
external information on the CP phase pd- The most precise value of pd is determined 
from B^ —)■ J/ipK^ decays, but this determination is also affected by penguin pollution. 
A recent study of the penguin effects in B^ —>■ , B^ —)■ J/xjjK^ ^ B^ —)■ J/ipir^ and 

B^ —)■ J/ipKg decays is performed in Ref. j^, with the latest numerical update 
including the results from Refs. [6,55,^, leading to (pd = 0.767 ± 0.029rad. 

In addition, a second set of Hi observables can be constructed by replacing 
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^ J/p)K 


*0 


by B^ —)■ J/%pp^ in Eq. 
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To minimise the theoretical uncertainties 
associated with the use of these Hi observables, the strategy proposed in Ref. [I^ 
adopted. That is, the relation 


IS 




^ J/i. 4 ,) 


A'M 


Ai 


A(R0 ^ j/PjK*^) 


A(R0 



(37) 


between the hadronic amplitudes in —>■ J/ipK*^ and B^ —)■ J/pjp^ is assumed, and 

therefore relying on theoretical input from LCSR is no longer needed. Instead, the ratio 
can be determined directly from the fit, providing experimental information on this 
quantity. Effectively, the three CP asymmetry parameters Ci and Si determine the 
penguin parameters and 9i. Thus, this result for a* and 9i predicts the values of the two 
observables Hi{B^ —)■ J/%pK*^) and Hi{B^ —)■ J/%pp^). By comparing these two quantities 
with the branching fraction and polarisation information on B^ —>■ J/'ipK*^, B^ —>■ J/p) p^ 
and —)■ J/ip<p, the hadronic amplitude ratios |A.'/A.j| can be determined. The impact of 

the Hi observables on the penguin parameters a* and 9i is negligible in the combined fit. 

For the combined analysis of —?• J/pj K*^ and B^ —)■ J/pj a modified least-squares 

fit is performed. External inputs on 7 = (73.2PPq)° and pd = 0.767 ± 0.029rad 
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are included as Gaussian constraints in the fit. The values obtained from the £t are 


ao = , 

«|| = 0.07^°;^^ , 
ax = 0.04+°;i2 , 


^0 = - {SStllsY ^ 
e\\ = -{85tliy , 
e^= (381-)°, 


.4,0 

^jl 

A\ 

A 


1 iqci+0.074 

1.100_QQ5g , 

1 9QQ+0'104 
r-^oo_Q QgQ , 


1 042+°-°®^ 


with the two-dimensional conhdence level contours given in Fig. which also shows the 
constraints on the penguin parameters derived from the individual observables entering 
the fif different bands. Note that the plotted contours for the two H observables do 
not include the uncertainty due to |4,V4.|. 

The results on the penguin phase shift derived from the above results on a* and 9i are 


(stat) (syst) rad , 

= 0 . 00 llQ;gj ;4 (stat) ±0.008 (syst) rad , 
^ 0.003lloj;4 (stat) ±0.008 (syst) rad . 


These results are dominated by the input from the CP asymmetries in —)■ and 

show that the penguin pollution in the determination of 0* is small. 


10 Conclusions 


Using the full LHCb Run I data sample, the branching fraction, the polarisation fractions 
and the direct CP violation parameters in J/l> K*^ decays have been measured. The 

results are 


B{B° J/'4)K*Y = (4.14 ± 0.18(stat) ± 0.26(syst) ± 0.24(/rf//,)) x 10"^ 


/o 

_/|| 


0.497 ± 0.025 (stat) 
0.179 ± 0.027 (stat) 
-0.048 ± 0.057 (stat) 
0.171 ± 0.152 (stat) 
-0.049 ± 0.096 (stat) 


± 0.025 (syst) 
± 0.013 (syst) 
± 0.020 (syst) 
± 0.028 (syst) 
± 0.025 (syst) , 


which supersede those of Ref. [1^, with precision improved by a factor of 2 — 3. The shift 
due to penguin pollution is estimated from a combination with the B^ —)• J/ipp^ 
and is found be to compatible with the result from the earlier analysis. 
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Appendices 


A Angular acceptance 


To take into account angular acceptance effects, ten normalisation weights, are com¬ 


puted and embedded in the normalization integral of the angular distribution given in 
Eq. 1^ following the procedure described in Ref. 57 
basis, the htting PDF can be written as 


Using the transversity amplitude 


dr(6'ic, 6^^, (ph) 

dD 




X]z<fc I H^'h) ^'ay V'h) 


(38) 


where the real or imaginary angular functions FijiOxi 0^, p>h) are obtained when combin¬ 
ing Eq. 1^ and Eq. |5]j^ and where eQ{9K,0fj,,(ph) denotes the angular acceptance. The 
normalization weights correspond to the integrals 


^ _ I / 'Tl^[Fij{0K, dfi, ^h)] 9fj_, ifh) dD, if Fij G M , 

ifh)] eniOK, 9^, cpn) dD, if F,,- G I. 

In the absence of acceptance effects, the normalisation weights related to the interference 
terms are equal to zero by definition, whereas those related to each polarisation amplitude 
squared are equal to unity. Eight sets of normalisation weights are calculated separately, 
one for each 'mK-n+ bin and kaon charge. 

In order to correct both for imperfections in the detector simulation and for the 
absence of any S-wave component in the simulation sample, the weights are refined using 
an iterative procedure where the angular acceptance is re-evaluated recursively until it 
does not change signihcantly. Table gives one set of normalisation weights after the 


Table 6: Corrected angular acceptance weights for K vr"’' events lying in the first mx-^+ bin. 
The ^ij weights are normalised with respect to the ^oo weight. 



ij 

o 

o 

1 

(00) 

1.000 

2 

(HID 

+1.379 ±0.029 

3 


±1.388 ±0.030 

4 

(11^) 

±0.035 ±0.019 

5 

(Oil) 

-0.003 ±0.012 

6 

(01) 

±0.010 ±0.011 

7 

(SS) 

±1.190 ±0.019 

8 

(S||) 

-0.042 ±0.017 

9 

(ST) 

±0.029 ±0.016 

10 

(SO) 

-0.929 ±0.024 
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iterative procedure. The effect of this correction is below one standard deviation for all 
the normalisation weights except for the (SO) weight. This is expected due to the rapid 
efficiency drop close to cos 9 k = 1 which directly impacts the (SO) weight. At each step 
of this procedure the simulation sample is corrected both for the absence of an S-wave 
component and for the imperfections in the detector simulation. For the first correction, 
the angular fit result to data is used, whereas for the second the kaon and muon track 
momentum distributions of data are used. In both cases the correction is implemented by 
assigning weights to each event of the simulation sample. 


B Correlation matrix 


The statistical-only correlation matrix of the angular parameters obtained from the fit to 

is given in Table Here, the superscript / = 0,1,2, 3 in Fg 


data, as described in Sect. 8.1 


and (5g represent the number of the mK-n+ bin as defined in Table 
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Table 7: Statistical correlation matrix for the parameters from the angular fit 
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